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ABSTRACT 



We analyze a 96 ks Chandra hetgs observation of the single G-type giant HR 9024. 
The high flux allows us to examine spectral line and continuum diagnostics at high temporal 
resolution, to derive plasma parameters. A time-dependent ID hydrodynamic model of a 
loop with half-length L = 5 x 10 11 cm (~ R+/2), cross-section radius r = 4.3 x 10 10 cm, with 
a heat pulse of 15 ks and 2 x 10 11 erg cm" 2 s" 1 deposited at the loop footpoints, satisfactorily 
reproduces the observed evolution of temperature and emission measure, derived from the 
analysis of the strong continuum emission. For the first time we can compare predictions 
from the hydrodynamic model with single spectral features, other than with global spectral 
properties. We find that the model closely matches the observed line emission, especially for 
the hot (~ 10 s K) plasma emission of the Fexxv complex at ~ 1.85A. 

The model loop has L/R+ ~ 1/2 and aspect ratio r/L ~ 0.1 as typically derived for 
flares observed in active stellar coronae, suggesting that the underlying physics is the same 
for these very dynamic and extreme phenomena in stellar coronae independently on stellar 
parameters and evolutionary stage. 



1. Introduction 

The first X-ray stellar surveys sho wed widespread pre sence of coronal emission in the 



cool half of the H-R diagram (see e.g., IVaiana et al.lll98ll for a review). Improved spatial 



and spectral resolution have provided us with more powerful tools for investigating the char- 
acteristics of the X-ray coronal activity in late-type stars, and for exploring the underlying 
processes at work in stars with different stellar parameters, and evolutionary stage. 



1 Massachusetts Institute of Technology, Kavli Institute for Astrophysics and Space Research, 70 Vassar 
street, Cambridge, MA 02139, USA; testa@space.mit.edu 

2 Dipartimcnto di Scicnzc Fisichc & Astronomiche, Sezione di Astronomia, Universita di Palermo Piazza 
del Parlamento 1, 90134 Palermo, Italy 

3 INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 

4 Imperial College London, Blackett Laboratory. Prince Consort Rd, London, SW7 2AZ. UK 



-2- 



X-ray observations of late-type stars show that coronal phenomena observed at close 
range on the Sun are common in late-type stars though they occur on more extreme scales 
in very active stars. For insta nce, stars at higher activity levels have reach much higher tem- 



peratures and densities (e.g., ISanz-Forcada et al.ll2002l ; iTesta et al.l 12004 iNess et all 12004 ) 



than typically observed for solar coronal plasma. Analogously, coronal flares are routinely 
observed in late-type stars with characterist ics similar to those observed in solar flares, e.g., 
fast rise and slow decay (e.g., iReald 120021 for a review), but at the same time the flare 
frequency and intensity can be dramatically larger than observed on the Sun, with stellar 
X-ray luminosity i ncreasing orders o f magnitude with respect to the quiescent level (e.g., 
Favata et allboool : bsten et alibood ). 



The analysis of lightcurves during flares provides us with insigh ts into the characteristics 

of the magnetic field (e.g 



Schmitt fc Favata 1999 



Favata et al. 


2000; 


Reale et al. 


2004) 



great deal of information on the coronal heating and on the plasma structure morphology can 
be inferred from the detailed modeling of stellar flares; for instance, if enough data statistics is 
available for moderately time-resolved spectral analysis, the study of the complete evolution 
of a flare allows to infer whether the flare occurs in closed coronal structures (loops), what 
is the size of these flaring structures, whether continuous heating is present thr oughout the 



flare, and even constraints on the location and distribution of the heating (see iReale et al. 
2004 and §4.21 for a detailed discussion). 



This work presents a detailed modeling of a large X-ray flare on the single evolved 
Gl III giant HR 9024, observed with the Chandra High Energy Transmission Grating. Most 
of the flare is observed, from the rise phase to the late decay, allowing to constrain a detailed 
hydrodynamic simulation of the flaring structure. The analysis of this flare is especially 
interesting in the context of X-ray activity of evolved giants. HR 9024 is an intermediate 
mass star (M* ~ 2.9M ) in the Hertzsprung gap, that is in its initial rapid (< 1 Myr) phase 
of post-main sequence evolution when it enters the c ool region of the H-R diagram, and it 
develops a subphotospheric convective layer (see e.g.. iPizzolato et al.ll2000l ). The intermedi- 
ate mass Hertzsprung gap giants are strong X-ray emitters (Lx ~ 10 31 erg s" 1 ) while their 
main-sequence progenitors, late-B or early-A type, are X-ray dark lacking the fundamental 
ingredients to sustain X-ray activity, either magnetic-dynamo as in late-type stars or strong 
winds as in massive stars. The onset of an efficient convective layer, together with their 
typically fast rotation rates (due to the little if any loss of angular momentum in their main 
sequence phase) is thought to generate a dynamo mechanism sustaining the X-ray activity of 
these evolved yellow giants . The young coronae of these stars a re characterized by high tem- 



perature and density (e.g.. lAyres et al 



1998 



stars at the same high activity levels (e.g. 



Testa et al 



20041) . similarly to low-mass a ctive 



Sanz-Forcada et all 12002 ; ITesta et al.l 120041 ); on 
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the oth er hand, there is ev idence for significant differences such as much lower coronal filling 
factors fjTesta et al.ll2004l ). and very limited flaring activity. The evolved intermediate mass 
giants, both in the Hertzsprung gap and in the p ost helium flash " clump" (the relatively 
long-lived, ~ 70 Myr, core Heliu m burning phase 



stant X-ray emission level (e.g., lHaisch &: Schmitt 



Audard et al. 



2004 



Scelsi et al. 



sources ( Haisch fc Schmitt 



1994 



Avres et alJllgggL sh ow extremely con- 



1994 



Gondoin 



2003; Testa et al. 



2004 



20041) and only a few flares have been observed on these 



Ayres et al.lll999l . and this work). 



The opportunity of deriving a loop length from the analysis of the flare is interesting also 
in that it gives us an opportunity to probe the structuring of the corona, which in principle 
can be significantly different from dwarf stars coronae. In fact in these giants the gravity 
is considerably lower than for MS late-type stars therefore yielding a larger scale height, 
and possibly allow ing the existence of very extended coronae (as suggested for example by 
Ayres et al.ll2003l ). In the case of HR 9024 the surface gravity is only ~ 0.02g Q implying a 
scale height of about 3 stellar radii at 10 7 K (30-R* at 10 8 K). 



HR 9024 (HD 223460, OU And) is a moderately rotating (v sinz ~ 21 km s^. lde Medeiros fc Lebre 



19921 ) chromospherically active si ngle giant, not too wel l studied even though it is a bright 
and close by object (d ~ 135 pc, iPerryman et al.lll997l ). The stellar parameters are listed 
in Table [TJ Several X-ray observ ations of HR 9024 e xist indicating high and constant X-ray 
luminosity of a few 10 31 erg s _1 . ISingh et al.l (119961 ) analyzed ROSAT PSPC observations of 
chromospherically active stars and derived an X-ray luminosity 4.3 x 10 31 erg s _1 in the 
0.2 — 2.4 keV energy band; this transla tes to Lx~ 2 . 6 x 10 31 erg s _1 using the revised distance 
of 135 pc instead of 175 pc assume d by ISingh et al.l (119961 ). This value is in a good agreement 
with values obtained by iGondoinl (120031 ). in the 0.3 — 2 keV range, from two short XMM- 
Newton observations of HR 9024 showing little if any variability: 2.7 x 10 31 erg s" 1 

(Rev. 107, 6 ks exposure time), and Lx~ 2.3 x 10 31 erg s" 1 (Rev. 200, 3 ks exposure time); the 



total X-ray luminosities in the 0.3- 
respectively. 



10 keV range are 3.8 x 10 erg s and 3.0 x 10 erg s 



Table 1: Stellar parameters. 



Spec, type 


d 


R 


M 


log L bol P rot 




[PC] 


[Re] 




[erg s^ 1 ] [days] 


Gl ill 


135 a 


13.6 b 


2.9 C 


35.4 d 23.25 c 



Refere nces: (a) from SIMBAD; (Vl lSmgh et aD <jl99fih : (c) IGondoinl (Il999h : fd nFlowerl lll99fih : (e) ISingh et al 
dl99fih . 



Our approach here is to inspect the flare lightcurve and spectra and derive some quan- 
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tities relevant to set up a loop hydrodynamic model and to constrain the initial parameters. 
The numerical solution of the hydrodynamic equations allow us to synthesize in detail the 
emission (the so-called forward modeling) as it would be really observed and therefore to 
compare directly to the data. This allows us to have a feedback on the model and to re- 
fine the model parameters. Once obtained a good description of the global features, i.e. 
the lightcurve and the overall evolution of the temperature and total emission measure, the 
model will have constrained the loop length and aspect, and the heating function. The com- 
parison with the fine details of the data analysis will give us insight into the flare density and 
thermal structure and evolution, and help us to interpret the results of the data analysis. 

The paper is structured as follows: we describe the observations in £J2, and the methods 
both for the spectral analysis and the hydrodynamic modeling in §121 in §H we present the 
results of our study, and in §|5] we draw our conclusions. 



Observations 



We ana lyzed the Chandra High Energy Transmission Grating Spectrometer (see lCanizares et al. 



20001 . 120051 for a description of the instrumentation) observation of the X-ray active single 
Gl giant HR 9024. The data were obtained from the Chandra Data Archive^ and have 
been reprocessed using standard CIAO 3.3 tools and analysis threads. Effective areas and 
line responses (ARFs and RMFs) were calculated using standard CIAO procedure^. The 
characteristics of the hetg observation are listed in Table [2j 



Table 2: Parameters of hetg observation. 



ObsID 


Start 


End 


^exp 

[ks] 


logL X a log L X , pea fc b 

[erg s" 1 ] [erg s" 1 ] 


1892 


2001-08-11 
00:19:04 


2001-08-12 
03:43:14 


95.7 


31.8 32.1 



a From MEG spectrum, in the 2-24A range. Averaged over the whole observation. 

"Peak X-ray luminosity obtained from the MEG spectrum integrated in the time interval 10-15 ks from the 
beginning of the observation. 

Figure [T] shows the Chandra High Energy (heg) and Medium Energy (meg) Grat- 
ings spectra for the ~ 96 ks observation of HR 9024. The spectrum is characterized by a 



5 http://cxc. harvard.edu/cda 



http: / / cxc.harvard.edu/ ciao/threads/ 
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strong continuum emission, indicating high temperature plasma, as well as the unusually 
strong Fexxv complex at ~ 1.85A (logT[K] ~ 7.8), and Caxx resonance line at 3.02A 
(logT[K] ~ 7.7); these characteristics are some what extreme wh en compared with the spec- 
tra of other very active stellar coronae (see e.g. iTesta et al.ll2004l ). Other prominent features 
in the spectrum are the H-like and He-like lines of Si (~ 6.2, 6.7A respectively) and Mg 
(~ 8.4, 9.2A), and the strong emission lines of highly ionized Fe around llA. 
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Fig. 1. — Chandra heg (red) and meg (blue) spectra obtained in a 96 ks observation of the 
single giant HR 9024. Line identification for many prominent spectral features is provided. 



The lightcurve of the summed heg+meg dispersed photons integrated in the 1.5-26A 
range, presented in Figure [H shows clear variability: the X-ray emission level rises steeply 
at the beginning of the observation by a factor ~ 3 in about 15 ks; after a slow decay 
on a timescale of about 40 ks the lightcurve rises again, with peak around 80 ks from the 
beginning of the observation. The lightcurves for a hard (1.5-12A) and a soft (12-26A) 
spectral band, also shown in Fig. [2J indicate hardening of the spectrum corresponding to the 
two peaks of the emission, typical of stellar flares. The peak luminosity above 10 32 erg s _1 
is extremely high when compared with typical energies of stellar flares observed in active 
coronae of late-type stars. 
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2.0 




time [ks] 

Fig. 2. — Lightcurve obtained as the sum of total counts of HEG and meg dispersed spectra, 
using a temporal binsize of 2 ks. Lightcurves in hard (blue) and soft (red) spectral bands are 
shown. The vertical green dashed lines mark the time intervals selected for the temporally 
resolved spectral analysis and the hydrodynamic modeling (see $3]). 



3. Analysis 

Spectra were analyzed with the PINTofALE0 IDlfj software (IKashyap fc Drake! l2000l ). 
The high resolution spectra provide several plasma diagnostics (temperature, density, abun- 
dances, emission measure distribution), from the analysis of both continuum and emission 
lines, and from the lightcurves in different spectral bands or in single lines. The high flux 
allows us to examine spectral line and continuum diagnostics at high temporal resolution. 



3.1. Spectral analysis 

Continuum emission: The strong continuum emission in this spectrum provides us with 
constraints on the hot plasma on which it strongly depends. Specifically, temperature and 
emission measure, EM, can be estimated through a fit of the continuum, and their evolution, 
probing the hottest plasma component, allows us to constrain the hydrodynamic model of 
the flaring structures (see §3.21 and §4. 2ft . In order to derive T and EM, we fit the continuum 



7 http: / /hea- www.harvard.edu/PINTofALE 
interactive Data Language, Research Systems Inc. 
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simultaneously in the heg and meg spectrum, selecting narrow spectral regions which can 
be assumed as "li ne-free" to a reaso nable extent, on t he basis of predictions of the atomi c 



databases APED flSmith et al.l l200lh . and CHIANTI flDere et al 



19971 : lLandi et all 120061 ). 



We fit the continuum with an isothermal model computed in PINTofALE, using CHIANTI, 
and convolved with the hetg spectral response. This model contains all the contributions 
to the continuum (free-free, free-bound, 2 photon), however at high temperature (logT[K] > 
7.5), as observed in this spectrum, the bremsstrahlung continuum is by far the dominant 
process for the formation of the continuum. The "line-free" spectral regions used for the 
fit are: 2.00 - 2.95A, 4.4 - 4.6A, 5.3 - 6.0A, 7.5 - 7.8A, 12.5 - 12.7A, 19.1 - 20A. The fit 
also provides an estimate for EM from the normalization parameter. The analysis of the 
strong continuum emission allows us to derive the evolution of the temperature and emission 
measure during the flare with high temporal resolution. The time intervals chosen for this 
analysis are shown in Fig. [2] superimposed on the lightcurve; they have been selected in order 
to have high temporal resolution but also good constraints on the plasma parameters. 



Line em ission: Th e line fluxes are determined using the technique of spectral fitting de- 



scribed in lTesta et al.l (120041 ) . The measured line fluxes have been used for the reconstruction 
of the emission measure distribution and the determination of abundances ( §4.ip . and for 
the comparison with the results of the hydrodynamic model ( §4.2p . 



Abundances and Emission Measure Distribution: The emission measure distribu- 
tion, EM(T), is derive d through a Markov-Cha in Monte-Carlo analysis using the Metropolis 

algori thm (MCMC[M]: lKashyap fc Drakdll998l ) on a set of line flux ratios, as in lGarcia- Alvarez et al 
(120061 ). O lines are the coolest lines used and Ar lines the hottest, covering a temperature 
range logT[K] ~ 6.2 — 7.8. Coronal abundances are evaluated on the basis of the derived 
EM( T): the abundance is a sca ling factor in the line flux equation to match the measured 
flux (jGarcia- Alvarez et al.ll2006l ). 



3.2. Hydrodynamic modeling 

The inspection of the light curve and of the evolution of the temperature and of the 
integrated emission mea sure (EM) during the flare allows to set up a detailed model of 
the flaring structure(s) (IReale et al.l I2004J ). as described in detail in §4.21 We derive the 
temperature and EM parameters through the analysis of the continuum emission as described 
above ( §3.ip . the continuum emission being strong and probing the hottest plasma. 
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The ID hy drodynamic mode l , solves time-dep endent plasma equations with detailed 
energy balance fjPeres et al.l Il982l ; iBetta et al.l 119971 ) , with a time- dependent heating; func- 
tion defining the energy release triggering the flare (see e.g., iReale et al.lll997l . 120041 ) . The 
coronal plasma is confined in a closed loop structure, where plasma motion and energy trans- 
port occur only along magnetic field lines. For the i nitial atmosphere we assume a loop in 
hydrostatic equilibrium and detailed energy balance (ISerio et al.lll98ll ) with maximum tem- 
perature (at the apex) T max = 2 x 10 7 K; the initial conditions do not influence much the 
evolution of the plasma after a very short time. 

The outputs of the hydrodynamic simulations are distributions of temperature and 
density along the loop sampled at regular times throughout the flare evolution. From these 
plasma parameters we synthesize the corresponding hetg spectrum of each (isothermal) 
plasma volume along the loop at each given time, using isothermal models folded with 
the hetg spectral response. We then integrate along the loop to obtain the overall hetg 
spectrum of the multi-thermal plasma in the flaring structure. 



4. Results 

4.1. Emission measure distribution and abundances 

As a first step we carry out an analysis of the spectrum to study the global characteristics 
of thermal distribution and abundances. To obtain this information we need strong lines, 
from different elements, covering a wide temperature range. We have to integrate on time 
intervals long enough to have good photon statistics for each line. This limits our temporal 
resolution and we are not able to derive EM(T) and abundances in several portions of the 
flare, but only in two different portions of the observation: during the flare (i.e. using the 
spectrum integrated over the first 40 ks of the observation), and outside the flare (40-96 ks). 
We note that hereafter we will label the parameters derived outside the flare as " quiescent" ; 
even though the corona seems to undergo another dynamic event this second flare is on a 
smaller scale with respect to the first one. The analysis of the abundances is useful also 
for the synthesis of the spectra from the results of hydrodynamic modeling and a consistent 
comparison with the data. 

Table [3] lists the fluxes of the spectral lines used for this analysis, measured in the two 
phases of flare and quiescence. 

Fig. [3] and H] show the emission measure distribution and the coronal abundances as 
derived from the flare spectrum, and from the quiescent emission. 
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Table 3. Measured fluxes (in 10 6 photons cm 2 sec 1 ) of the spectral lines used in our 
analysis, during the flare and outside the flare. 



ion 


Aobs |AJ 


log i max [KJ 


IlUXflarc 


tlux 


quiesc 


use 


Fexxv 


1.853 


7.8 


94 ± 19 


13.5 


± 


8.4 


Al 


FeXXV 


1.864 


7.8 


49 ± 16 


16.0 


± 


8.4 


AI 


Caxx 


3.024 


7.7 


21.7 ± 6.8 




< 


4.4 


s 


Caxix 


3.187 


7.4 


18.7 ± 6.4 


5.8 


± 


3.7 


s 


Ar xvni 


3.734 


7.6 


14.1 ± 7.0 


4.2 


± 


3.6 


s 


Ar xvn 


3.945 


7.3 


13.5 ± 6.7 


4.7 


± 


4.0 


s 


Sxvi 


4.729 


7.4 


20.3 ± 8.7 


10.8 


± 


5.7 


s 


Sxv 


5.041 


7.2 


30 ± 11 


9.5 


± 


6.0 


s 


Si xiv 


6.183 


7.2 


82.9 ± 6.0 


42.2 


± 


3.6 


SM 


Si xni 


6.648 


7.0 


33.1 ± 4.2 


25.7 


± 


3.0 


s 


Al xni 


7.170 


7.1 


15.9 ± 4.5 


8.1 


± 


3.1 


s 


Al xn 


7.759 


6.9 


6.2 ± 3.8 


5.6 


± 


3.1 


s 


Mgxn 


8.422 


7.0 


80.7 ± 6.4 


47.5 


± 


4.5 


SM 


Mgxi 


9.168 


6.8 


31.0 ± 4.7 


15.2 


± 


2.8 


s 


Ncx 


12.132 


6.8 


142. ± 12. 


113. 


± 


9.3 


s 


Fexxi 


12.284 


7.0 


79. ± 11. 


53.5 


± 


7.2 


s 


Ne ix 


13.448 


6.6 


20. ± 10. 


22.6 


± 


7.4 


s 


Fe xvni 


14.201 


6.9 


50. ± 15. 


30.1 


± 


9.1 


s 


Fe xvn 


15.012 


6.7 


107. ± 22. 


43. 


± 


13. 


s 


vin 


18.965 


6.5 


270. ± 42. 


179. 


± 


29. 


s 


Ovu 


21.602 


6.3 


51 ± 44. 




< 25 


s 


Nvn 


24.778 


6.3 


65. ± 50. 


72. 


± 40. 


s 



a Tempcraturc of maximum formation of the line. 

b Use indicates whether the line was used in the spectral analysis 
("S"), for the emission measure reconstruction and abundance deter- 
mination, and whether the feature was used for the direct comparison 
with the loop model ("M"; see $h2|). 
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The thermal distribution of the coronal plasma in HR 9024 appears dominated by hot 
(i.e., typical of flaring structures) plasma both during and outside the flare. Hot emission 
(T > 2 x 10 7 K) was found also from the analysis of two very short (6 ks, and 3 k s) XMM- 
Newt on observations of HR 9024 when the corona seems to be in its quiescent state (IGondoin 
20031 ). The main difference between the EM(T) of flare and quiescence is at the hot end 
of the temperature range, i.e. for logT[K] > 7.5, where the EM(T) of the flaring plasma is 
about one order of magnitude higher than the EM(T) of the plasma outside the flare. The 
EM(T) derived for HR 9024 presents interesting characteristics when compared to other 
active coronae. Specifically, the EM(T) of HR 9024 is characterized by a rather shallow 
slope, similar to the slope of hydrostatic loop models (EM(T) oc T 3//2 ; iRosner et al.lll978l ). 
whereas for several active coronae there is increasing evidence of steep, almost isothermal, 
emission measure distributions possi bly indicating the dynamic nature of the coronal loops 
composing them fjTesta et al.ll2005bl ): for instanc e, such steep EM ( T) are derived for other 



giants, e.g., the Hertzsprung gap giant 31 Com 



Griffiths &: Jordan 



Scelsi et all 120041). and the clump giants 



(3 Cet, dSanz-Forcada et al.ll2002l). and Capella (IDupree et al.lll993f). and for o ther active 
stars (e.g., 



19981 : brake et al.ll2000l : ISanz-Forcada et all 120021 ). 



E 

o 

o 
in 
O 




6.5 7.0 7.5 
log (T) 




6.5 7.0 7.5 
log (T) 



Fig. 3. — Emission measure distribution derived from the flare portion of the spectrum (0- 
40 ks; left), and from the spectrum outside the flare (40-96 ks; right). The comparison shows 
that the main difference between flare and quiescent emission resides in the hot end of the 
temperature range, i.e. for logT[K] > 7.5 where the EM(T) of the flaring plasma is about 
one order of magnitude higher than the EM(T) of the plasma outside the flare. 



The coronal abundance pattern observed in HR 9024 (Fig. HJ) pr esents characteristics 



simil ar to other intermediate mass giants, i.e. little or no FIP effe ct (IGarcfa-Alvarez et al. 



20061 ). All abundances are plotted relative to the solar mixture of lAsplund et al.l (120051 ) as 
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photospheric abundances for HR 9024 have not been determined. We find significant changes 
in the coronal abundances during the flare with respect to the quiescent phase: all elements 
appear to be enhanced, except possibly Ne, and the Fe abundance in particular is found to 
increase by almost a factor 3. 
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Fig. 4. — Coronal abundances relative to solar abundances (lAsplund et al.ll2005l ) as a func- 
tion of the First Ionization Potential (FIP) derived for HR 9024 during the flare (square 
symbols) and the quiescent (triangles) phase. 



Both abundances and EM(T) for HR 9024 have been derived also by iNordon Sz Behar 



( 12006I ). and their findings are in agreement with the results presented here. 



4.2. Flaring loop modeling 

High level diagno stics of the flare can be obtained from detailed hydrodynamic modeling 



of the flaring plasma (IReale et al.lll988l . I2004T ) . The analogy of stellar to solar flares suggests 



that stellar as solar flares mainly occur in closed magnetic structures (coronal loops). It is 
customary to assume that the bulk of the flare involves a single coronal loop, which can be 
then investigated through the analysis of the flare characteristics. 

The plasma confined in coronal loops can be described as a compressible fluid which 
moves and transports energy exclusively along the magnetic field lines. We can then model a 
flare by solving the time-dependent hydrodynamic equations in the coordinate along the loop 
with the assumption of an appropriate transient input heating function. The main model 
parameters include the loop length and the intensity, distribution and duration of the heat 
pulse. These parameters are constrained by the observed flare evolution, in particular in the 
decay phase, and by parameters derived from the data analysis, such as the temperature and 
emission measure at flare maximum and the timing of the maximum. Scaling laws such as 
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those shown in iReale et al.l (119971 ) help us to constrain the loop length from the light curve 
decay time and the flare maximum temperature: 



TxVT7 



^ 120/(C) (1) 
where Lg is the loop half length (10 9 cm), tx the decay time of the light curve, T-j the 
flare temperature (10 7 K), /(£) > 1 a correction factor which takes into account possible 
significant heating during the decay, that might make the decay slower. 

Other parameters, such as the loop initial atmosphere, or, for low gravity stars, the stel- 
lar surface gravity, have much less influence on the flare evolution. We compute the evolution 
of the loop plasma by solving the time- dependent hydrodynamic equation of m ass, momen- 



tum and energy conser vation for a compressible plasma confined in the loop (IPeres et al. 



19821 ; iBetta et al.l 119971 ). including the relevant physical effects such as the plasma ther- 
mal conduction and radiative losses. The gravity component along the loop is computed 
assuming a radius i?* = 3R Q and a surface gravity g* = 0.1g Q . 

We integrate numerically the equations over a time range of 50 ks. From the density and 
temperature distribution of the plasma along the model loop, we synthesize the expected 
plasma X-ray model spectrum, and finally convolve with the Chandra/ 'meg spectral response 
to produce synthetic counts spectra as described in §3.21 



The flare light curve integrated in the whole Chandra/ meg band resembles qu ite closely 
the X-ray light curves of flares recently observed on PMS stars ( iFavata et al.ll2005l ). For this 
reason, we have assumed as a first set of parameters the same values which best describe a 
specific flare observed in the Orion region, i.e. a loop with constant cross-section and half- 
length L = 10 12 cm, symmetric around the loop apex. The flare simulation is triggered by 
injecting a heat pulse in the loop which is initially at a temperature of ~ 20 MK. Two heat 
pulses are deposited with a Gaussian spatial distribution of intensity 10 erg cm -3 s -1 and 
width 10 10 cm (1/100 of t he loop half-length ) at a distance of 2 x 10 10 cm from the footpoints, 
i.e. very close to them (IReale et al.l 120041 ) . After 20 ks the heat pulses are switched off 
completely. With the above parameters, the light curve is reproduced with good accuracy 
but we have noticed that the temperature values obtained from fitting the spectra with 
isothermal plasma emission m odels are significant ly higher than those derived from the data 
all along the flare evolution (ITesta et al.l l2005al ). Therefore, we have refined the model 
parameters to improve the fitting: we have reduced the input heating enough to account for 
the desired temperature decrease. From Eq.([T]), in order to maintain the same decay time, a 
temperature reduction implies a shorter loop length. After this feedback, we have obtained 
best results with a loop half-length L = 5 x 10 11 cm (half as before) and a heat pulse of 
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4 erg cm 3 s 1 (peak volumetric heating) lasting 15 ks (5 ks less than before). The total 
energy flux rate is 2.0 x 10 11 erg cm -2 s _1 . 
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Fig. 5. — Evolution of the temperature, density and pressure at the loop apex and of the 
maximum plasma velocity according to the loop hydrodynamic model of the flare. 



The computed evolution largely resembles the evolution computed for other stellar flares, 
although on larg er scales than for typical flares observed in late-type stars coronae (e.g. 
Reale et al.ll2004l ). Fig. [5] shows the evolution of the temperature, density, pressure at the 
apex of the loop and of the maximum plasma velocity. The heat pulses make the loop plasma 
heat rapidly (~ 3 ks) at temperatures about 150 MK and expand dynamically upward from 
the chromosphere at speeds about 1800 km s _1 to reach the loop apex also in about 3 ks. 
After this first impulsive phase, the temperature slowly decreases to about 130 MK while 
the heat pulse is on, and the evaporation of plasma from the chromosphere - which brings 
the density increase shown in Fig. [5] - continues substantially but less dynamically, with 
plasma speeds below 500 km s" 1 after 9 ks, when the plasma pressure becomes higher than 
1000 dyne cm -2 . The plasma pressure and density reach their peak of 1500 dyne cm -2 and 
4.5 x 10 10 cm -3 slightly (about 3 ks) later than the end of the heat pulse. Then they begin 
to decrease gradually with significant quasi-periodic fluctuations, while the temperature 
initially drops and then decreases more gradually. 

Figure [6] shows the integrated light curve assuming a loop cross- section radius of 4.3 x 
10 10 cm, i.e. 8% of the loop half-length. The spectra obtained from the loop modeling have 
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also been fit with isothermal models. The resulting evolution of the emission measure and of 
the temperature, and the T vs. EM diagram are shown in the figure. A visual comparison of 
the model results with the data and the data fitting results indicates that the loop simulation 
results are in good agreement with the data. 



1.2 




100 



20 40 60 80 100 27.10 27.20 27.30 27.40 27.50 

t [ks] 0.5 log(EM [cnT 3 ]) 

Fig. 6. — Comparison of observed lightcurve (top left), T (top right), EM (bottom left), 
and T-n (bottom right) evolution, with the corresponding quantities synthesized from the 
hydrodynamic loop model (solid thick lines) characterized by half-length L = 5 x 10 11 cm, 
and footpoint impulsive heating. Evolution of temperature and emission measure (top right, 
and bottom left panels) are derived from the fit with isothermal models to the continuum 
emission in the spectra integrated in the time intervals shown in Fig. [2j 

The dynamic model described above closely matches the evolution of temperature and 
emission measure derived from the continuum emission, and the lightcurve integrated over 
the entire wavelength range. The high quality data allows us to test the model and perform a 
detailed comparison of its characteristics with the observations. In Figure [7| we compare the 
meg spectra synthesized from the model with the observed spectra in 4 different phases of 
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the flare: rise (0-10 ks from beginning of observation) , peak (10-15 ks), early decay (15-30 ks), 
and late decay (30-45 ks). Figure [7] shows a general good agreement of the model spectra with 
the observations; however, some systematic discrepancies are observed. For example, in the 
spectrum of the late decay phase several lines, such as the Sixiv Lya (~ 6.2A), the hot Fe 
lines around llA, and the Ovill Lya (~ 19A), are overpredicted. There are several factors 
that may determine these departures from the observations. For instance, the abundances 
used for synthesizing the model spectrum are the ones derived from the flare spectrum and 
are kept fixed. Since we observe significant changes in the abundances after the flare ( §4.11) 
they may change on short time scales and affect the emission in the spectral lines. Another 
process, not taken into account in our modeling, that might affect the emission in the lines on 
short time scale is non-equilibrium ionization; ho wever, this effect is e xpected to be relevant 



only at the very beginning of the flare (see e.g., lOrlando et al.ll2003l ). Furthermore, in our 
simulation we assume that all the observed emission comes from the flaring structure, i.e. 
the emission of the quiescent corona (or of other secondary flaring structures) is assumed to 
be negligible whereas it could in fact contribute significantly to the "cool" (logT[K] < 7.2) 
lines at least, as suggested by the derived emission measures during flare and quiescence 
(Fig. ED. 

The high quality of the data allows us to test the model to a much higher level of 
detail than before, for instance allowing a direct comparison with the lightcurves in single 
spectral features. The evolution in very hot spectral features emitted at the bulk plasma 
temperature of the flaring loop (i.e. logT[K] > 7.5) would represent more meaningful tests 
for the model. However, there are only a few very hot lines present in the spectrum (e.g., 
Fexxv and Caxx), and only the Fexxv complex has enough signal to provide temporal 
resolution. Beside Fe xxv, we selected two other emission lines with high enough flux, Si xiv 
and Mgxn Lya lines, and derived their lightcurves using the same time intervals used for 
the analysis of the continuum and marked in Fig. [2J Figure M compares the lightcurves 
in these spectral features derived from the observations with those synthesized from the 
hydrodynamic simulation, assuming a loop cross-section radius r = 4.8 x 10 10 cm. The line 
fluxes for Si XIV and Mg XII are obtained from the fit to the spectra integrated in each time 
interval. For the Fexxv complex the lower signal to noise ratio does not constrain the fit; 
however, since in the relevant wavelength range the continuum emission is rather small with 
respect to the line emission we use the lightcurve integrated in a small wavelength range 
(1.83-1.89A) using a 3 ks temporal bin. 

The model reproduces extremely well the observed temporal behaviour of the Fexxv 
line complex. The lightcurves of Sixiv and Mgxn derived from the observation and from 
the model lie close to each other up to the peak of the flare. After the peak the observed 
emission in the two lines experiences the dramatic drop predicted from the model, but on 
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Fig. 7. — Comparison of HR 9024 meg spectrum (black histogram) and the meg spectrum 
synthesized from the hydrodynamic model (red histogram) in four different phases of the 
flare, from top to bottom,: rise (0-10 ks from beginning of observation), peak (10-15 ks), 
early decay (15-30 ks), and late decay (30-45 ks). Spectra are shown in counts per second, 
and for better readability we split the spectral range in two plots showing the spectral range 
1.5-12A and 12-22A respectively. The small panel below each plot shows the x residuals. 



much shorter timescales, anticipating the model by slightly more than 10 ks. Other than the 
effects discussed above possibly causing this discrepancy, we can identify another effect that 
may play a role. In the model of the flaring loop the plasma contributing the most to these 
"cool" lines (logT[K] ~ 7.0 — 7.2) is located in the transition region of the loop, that is the 
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Fig. 7. — continued. 



region where the plasma temperature and density steeply change from the dense and cool 
chromosphere to the hot and rarefied coronal portion of the loop. In this region some models 
predict an opening of the magnetic flux tube which could have a substantiall y larger cross 



section in corona than in the low region connecting it to the photo sphere (e.g.. lGabriellll976 



Schrijver et al.lll989l ; iLitwin fc Rosnerlll993l ; ICiaravella et al. 



1996I ). In this possible scenario, 



the emission of the flaring structure in these cool lines would be substantially reduced and 
the resulting excess would be th erefore attributed to the background corona or to nearby 
secondary flaring structures (e.g.. lf3etta et al.ll200ll ). 
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Fig. 8. — Lightcurves derived from hetg spectra in selected spectral features (histograms 
with error bars) compared with the prediction from the hydrodynamic model (thick solid 
line), assuming a loop cross-section radius 4.8 x 10 10 cm. 

Finally, we compared the emission measure distribution of the loop model with the 
EM(T) derived from the flare spectrum. Figure [9] shows that also for the emission measure 
distribution the model predictions are in very good agreement with the observations. This 
agreement is a further indication of the validity of the single loop model for this flare and, on 
the other hand, of the EM inversion method. In our modeling of the X-ray emission during 
the flare we assume any persistent quiescent emission to be negligible. As mentioned in §4.1[ 
the Chandra observation does not provide a good determination of the quiescent emission as 
outside the large flare constituting the focus of our analysis the corona is undergoing another 
significant, though smaller, dynamic event (its average Lx ~ 4.2 x 10 31 erg s -1 is about twice 
the luminosity value of other X-ray observations; see §Q]and £J5]). Other X-ray observations 
of this corona suggest that the quiescent emission is in fact negligible to a large extent with 
respect to the large and hot flare observed with Chandra (see §5]). 



5. Discussion and conclusions 

Our analysis has shown that the observed evolution of the very hot coronal emission of 
the single giant HR 9024 is reproduced extremely well by a model characterized by loop semi- 
length L = 5 x 10 11 cm (~ R+/2), and impulsive footpoint heating triggering the flare. The 
heating pulse lasts 15 ks, and it is shifted by 8 ks preceding the beginning of observation. 
The peak volumetric heating is 4 erg cm" 3 s _1 , with a corresponding total heating rate 
po -^q33 er g g -i rj^g ^tailed analysis of the plasma evolution in the temperature-density 
diagram (Fig. [6]) put tight constraints on the flare morphology and on the characteristics 
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Fig. 9. — Comparison of emission measure distribution derived from the flare spectrum 
(histogram with error bars) with the EM(T) of the hydrodynamic loop model (solid thick 
line). The dashed line shows the EM(T) of the initial static loop with maximum temperature 
of 20 MK. 



of the heating. The delay of the EM peak with respect to the temperature peak indicates 
that the X-ray emitting plasma is confined in a closed structure, whereas in non-confined 
coronal flares the EM evol ves simultaneously with the temperature as shown by the analysis 
of hydrodynamic models (IReale et al.l I2002T ) . The steep slope of the decay path in the 
temperature-density diagram implies that after the initial phase when the energy is released, 
the loop undergoes a pure cooling evolution, and no addition heating is needed to explain 
the evolution of the flaring plasma. 

As discussed in §4.2! the initial phases of the flare are very dynamic and the plasma 
fills up the loop expanding upward at very high speed reaching ~ 1800 km s" 1 . After 
about ~ 10 ks, therefore corresponding to the start of the Chandra observation (considering 
the shift of 8 ks between the energy release and the start of the observation), the plasma 
speed rapidly decreases; however, during the first 5-10 ks of the observation plasma flows 
with speed up to a few 100 km s _1 are present. Such speed can be in principle resolved 
by hetg. The plasma characterized by such high speed is plasma at temperatures above 
60 MK, therefore, in order to investigate these speeds through Doppler shifts we have to 
search for shifts in very hot lines. However, such hot lines, e.g. Caxx, are at low A where 
the effective area is low and the instrument is not very sensitive. We searched for shifts in 
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the Caxx line but we did not detect any significant shift. We note that, even in the case 
the instrument sensitivity were not a limiting factor, the detection of line shifts is expected 
only for a preferential orientation of the loop with respect to the line of sight. 

The high photon statistics of the hetg data allowed us to derive light curves in a few 
relevant hot lines and, for the first time for a stella r flare, they could be compared in d etail 
to those predicted by a hydrodynamic model (see iPeres et al.l 119871 ; iBetta et al.l l200ll . for 
analogous comparisons for a solar flare). 

To a more general level, we could also compare the distributions of emission measure 
obtained from the data analysis and from the flare loop model and found a good agree- 
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From the normalization of the model lightcurves we derive an estimate of the loop cross- 
section and therefore of its aspect ratio a = r/L. We note that the loop cross-section radius, 
r, is a free parameter for which we can obtain independent estimates from the normalization 
of the different lightcurves: (a) integrated meg dispersed counts, (b) EM derived from the 
analysis of the continuum (Fig. [6]), and (c) fluxes of spectral features (Fig. [8]), therefore 
providing a cross-check of the consistency of our model. The obtained values of r = 4.9 x 
10 10 cm, 4.3 x 10 10 cm, and 4.8 x 10 10 cm respectively, all agree within 15%, and they imply 
a loop aspect ratio a ~ 0.09. Figure [TU] shows a sketch of the morphology of the flaring 
structure as inferred from our modeling of the observed Chandra-WETG spectra. 

As mentioned in §4.21 the loop model reproducing the observation has roughly the same 
parameters of models satisfyingly reproducing othe r large flares, specifically the flares ob- 
served in pre-main sequence stars (IFavata et al.ll2005l ). However, the derived loop parameters 
reveal some dissimilarity possibly indicating fundamental differences between the processes 
at work in the two cases. The flaring loop reproducing the flare on HR 9024 is characterized 
by loop semi-length L ~ R*/2, and aspect ratio r/L ~ 0.1, as typically found for low mass 
stars and on the Sun. In contrast, the modeling of the flare on the pre-main sequence stars, 
in particular the detailed model for COUP source 1343, yields a very large loop semi-length 
corresponding to s everal stellar radii, and with a much smaller aspect ratio r/L ~ 0.02 
(IFavata et al.ll2005l ). These findings suggest that HR 9024 is characterized by a "normal" 
corona, while in the young stars showing evidence of elongated X-ray emitting structures, 
possibly connecting the star to the disk, processes profoundly different may be at work. 

Large flares as the one observed in HR 9024 are very unusual in single intermediate 
mass giants, and there are only a few such events studied in detail. This study thus allowed 
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Fig. 10. — Sketch showing the geometry of the flaring loop of HR 9024 as derived from the 
hydrodynamic modeling of the flare: R* is the stellar radius, L is the loop semi-length, and 
r is the cross-section radius of the loop. 



us to investigate the flaring activity in a single giant, and provided us with a powerful tool 
to derive the size of its coronal structures. Even though the flaring structure modeled here 
is not necessarily representative of the typical coronal structures in X-ray active giants, 
as discussed above, its characteristics are consistent with the presence of a corona with 
properties typical of active coronae but scaled up to the larger stellar radius. The coronal 
properties derived through X-ray observations of these active giants might suggest a possible 
interpretation for the low frequency of flares. There is some evidence of very low coronal 
surfa ce filling; factors (< 10 -4 ) for these active giants as suggested from density analysis 



[e.g.. iTesta et al.ll2004j ). Considering the interaction of magnetic fields in active regions as a 
possible mechani sm increasing the X-ray flaring activity in active stars (IGiidel et al.lll997l ; 
Testa et al.ll2004l ). in these stars this kind of interaction might be much less frequent given 
to the very sparse presence of active region distributed on a large area. 

In this paper both for the modeling and its interpretation we have assumed that the 
quiescent emission is negligible with respect to the flaring emission. The Chandra observation 
does not provide a good grasp of the quiescent emission as the hardness ratio behaviour 
(and the emission measure distributi on) suggest s the presence of another dynamic event 



subsequent to the large flare (as in iReale et al.l I2004J ) . The XMM- Newton observations 



might represent the quiescent X-ray spectrum of HR 9024 not showing significant variability, 
however their very low exposure time (6 ks, 3 ks) allow us only to assess the variability on very 
short time scales. Our preliminary results of a ~ 50 ks new Suzaku observation of HR 9024 
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shows extremely constant X-ray emission and provides an estimate of the quiescent X-ray 
luminosity and temperature of about 2 x 10 31 erg s _1 and 3 x 10 7 K respectively (P. Testa et 
al. in prep.), consistent with the f indings based on XMM-Newton and ROSAT observations 
(IGondoinl 120031 ; ISingh et al.lll996h . This X-ray luminosity value which we assume being a 
good estimate for the quiescent emission of HR 9024, is about an order of magnitude lower 
than the peak luminosity of the flare observed with Chandra, therefore lending support to the 
assumption that the quiescent emission is negligible with respect to the flare emission. Also, 
the main results of our study rely on the emission of the very hot plasma which is completely 
dominated by the flaring structure, and therefore they should not be significantly affected 
by ignoring the much cooler quiescent emission. 

In this work we model the flare emission with a single loop structure and it is worth 
discussing whether and to what extent very different solutions can be ruled out. We model 
the bulk of the flare where presumably there is one dominant loop structure, as observed 
in many solar compact flares entirely occurring in single loops and even in more comple x 
flares, where one can often consider a dominant loop (e.g. lAschwanden &: Alexander! 120011 ). 
The single loop model satisfactorily reproduces the Chandra observations analyzed in this 
work and multiple loop models would not add insight and would include more free and 
unconstrained parameters. Furthermore, the single loop model is supported by the fact that 
: (1) the loop cross-section area is consistent with that of a single loop structure; (2) the 
thermal distribution is compatible with that of a single flaring loop; (3) the decay path in 
the density-te mperature diagram ha s quite a steep slope (~ 2) implying negligible heating 
in the decay fjSylwester et all 119931 ). in turn indicating a single loop structure , whereas 
flaring arcades are characterized by significant heating (e.g. iKopp &: Polettolll984l ). At later 
times of our observation (when the light curve rises again) other loops may be involved and 
become important in the evolution (e.g. iReale et al.l 120041 ) . The general expression of the 
loop length as a function of the observed decay time (Eq. [TJ also allows us to estimate the 
uncertainties on the derived loop length. In part icular, using the exp ression of the correction 
factor /(£) tuned for Chandra, as reported in iFavata et al.l (120051 ). we obtain Lg = 490. 
Considering typical uncertainties of the diagnostic formula, and the uncertainties on the 
observed temperature, we estimate an error of about 20-30% on the loop length, therefore 
yielding L = 4.9 ±1.5 x 10 11 cm. It is worth noting that this formula used for estimating 
the uncertainties on the derived loop geometry assumes uniform heating and provides only 
approximate values for the loop length. The hydrodynamic model provides us with a much 
higher level of detail, giving us diagnostics for the temporal and spatial distribution of the 
heating, for the loop aspect, and for the thermal distribution of the plasma. 
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